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JL We study Quantum Chromodynamics with eight flavours by use of lattice simula- 

tions and present evidence that the theory still breaks chiral symmetry in the zero 
temperature, continuum limit. This confirms that the lower end of the conformal 
window of QCD lies above N* = 8. 
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1 Introduction 



The pattern of chiral symmetry breaking in 577(3) gauge theories, and more 
specifically Q CD-like theories, with a large flavour content is not well un- 
derstood. Perturbative arguments suggest that chiral symmetry should be 
restored above a (yet to be determined) critical flavour value TVS, and below 
the value N\ F where asymptotic freedom disappears [lj. For flavour numbers 
between these values, the theory becomes conformal due to the presence of 
the Banks-Zaks infrared fixed point. 
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These findings were substantiated by a renormalization group estimate of the 
critical point [2J, and an upper bound on Nj has been proposed [3]. For 
a non-supersymmetric SU(N) theory with Nf massless flavours, one finds 



Nj < AN Jl — 16/81iV 2 . Instanton studies at large Nf [3] claimed a quali- 
tative change in behaviour at Nf = 6. Renormalization group flow equations 
predict Nf = 10(1) [3]. Recently, an all orders beta function inspired by the 
Novikov-Shifman-Vainshtein-Zakharov beta function of Af = 1 supersymmet- 
ric gauge theories has been computed, leading to a bound of the conformal 
window [6], in particular NJ > 8.25 for an SU(3) gauge group. These results 
have also been discussed from the more general perspective of higher dimen- 
sional representations [TJ. 

There is more than one reason for interest in these studies. First, the uncover- 
ing of the rich phase structure of vector- like theories such as QCD, is in itself 
a good step forward in our understanding of gauge field theories in general. 
Secondly, the extent and the nature of the conformal window is of interest 
to model builders, since it can be used to describe electroweak symmetry 
breaking. Finally, the exploration of the temperature dependence of the same 
phenomena adds one further perspective, as there should be a connection of 
some sort between the conformal phase of QCD and the quark gluon plasma 
phase. 

More recently, analytical studies of QCD at large Nf have been extended to 
finite temperature, and the critical line in the T, Nf plane has been predicted 
on the base of a truncated renormalization group flow calculation [8J. Note 
that the number of flavours can indeed be regarded as a continuous variable, 
since the Casimir's are polynomials in Nf j9]. This means that the critical line 
found in [8] is in fact a true phase boundary between a conventional Goldstone 
phase and a chirally symmetric phase, the zero temperature limit of which is 
the onset of the conformal phase. We sketch the resulting phase diagram in 
Fig- [D 

An interesting question is how the conformal phase of zero temperature, large 
Nf QCD is connected with the high temperature, small Nf quark gluon plasma 
phase. One important difference between the two is that the conformal phase 
possesses no inherent scale, while the other high T phase contains, of course, 
the temperature as a scale. As a consequence, the pion will be massless in 
the chiral limit of the conformal phase, while its mass will equal the lowest 
Matsubara frequency in the chiral limit at high T. Thus, the behaviour in the 
conformal window is in principle different from that in the high temperature, 
small Nf regime of QCD. However, it is unclear if the pion mass - or any 
other observable - can be regarded as an order parameter, or, in other words, 
if we should envisage a further critical line at the right side of the chiral line 
in Fig. [T] Additionally, the issue whether chiral symmetry restoration, for 
fermions in the fundamental representation, is always accompanied by, and 
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Fig. 1. Qualitative graph of the phase diagram in the T,Nf plane. 

coinciding with, a deconfinement transition, is an open and relevant problem. 
Much of the numerical evidence up to date points to their entanglement, but 
recent work advocates a separation of the two transitions, at least in the large 
N c limit jiopilj . It is very interesting to see what happens for large Nf. 

Lattice calculations provide a nonperturbative framework to attack this prob- 
lem. The phase diagram of Fig. [T] suggests a natural strategy for finding the 
lower limit of the conformal phase: simply, find for which Nf the transition 
disappears. As a by-product, the numerically determined critical temperature 
as a function of Nf can be compared with the analytic estimate of [8] , and the 
nature of the phase transition investigated. There is however one complication. 
Namely, it has been proposed - and to some extent proven - that the strong 
coupling limit of lattice QCD is always confining, regardless of the value of 
Nf (see again e.g. [H])- A simple argument is that, at large enough coupling, 
Nf means nothing. For instance, the statement that at Nf = 16 asymptotic 
freedom is lost, is a statement of the continuum weak coupling regime, which 
has no special relevance at strong coupling. Even for those Nf which are sup- 
posed to fall deeply in the conformal phase, a strong coupling limit still exists 
which breaks chiral symmetry and confines. 

This means that, even when the weak coupling continuum regime is conformal, 
a zero temperature 'lattice' phase transition will always take place between 
a strong and a weak coupling regime. In other words, for a fixed number of 



flavours, and for a given temporal extension of the lattice N t , we expect a 
phase transition to occur as a function of the lattice coupling to a 'strongly 
coupled' phase where chiral symmetry is broken, i.e. the chiral condensate 
is non zero in the chiral limit. This type of transition, called a bulk phase 
transition, has to be distinguished from a true thermal transition. 

We can thus envisage two types of lattice behaviour as a function of the 
number of flavours. When Nf is such that QCD in the continuum still breaks 
chiral symmetry and confines, i.e. when the continuum limit describes the 
ordinary phase of QCD, no transition is expected at zero temperature - or, 
equivalently, N t = oo, given that T = l/(N t a), with the lattice spacing a a 
function of the inverse lattice coupling (3 = 6/g 2 . If a lattice transition occurs it 
should get closer to (3 = oo and eventually disappear when N t — > oo. Instead, 
for the values of Nf for which an infrared fixed point exists, the lattice phase 
transition should survive the Nt —* oo limit, i.e. it is a bulk phase transition. 
The larger Nf, the closer this bulk phase transition should get to the infinite 
coupling j3 = limit. 

At a practical level, one should be able to distinguish between a bulk phase 
transition and a true thermal phase transition by studying the behaviour of the 
theory for a given Nf at different values of Nt and at varying lattice coupling 
(3. If the critical point remains non zero (3 C ^ when N t — > oo we have a 
bulk transition, thus providing a strong evidence that at the given Nf the 
continuum theory is chirally symmetric and falls into the conformal window. 
The alternative behaviour will point to a thermal transition happening at a 
given physical temperature T c . 

The authors of [12] performed a pioneering study at Nf = 10 and Nf = 
12 at small volumes (mainly 8 3 x 4) and relatively heavy masses, reporting 
interesting evidence for a persistence of a first order thermal phase transition 
at Nf = 10 and discussing its sensitivity to fermion masses. Subsequent work 
by the Columbia collaboration [13] suggested a rich structure of the phase 
diagram, dominated by lattice artifacts and was not completely conclusive. 
Iwasaki and collaborators performed a detailed lattice study, finding that a 
value as low as Nf = 6 should be set as the lower limit of the conformal 
window [14] . Note that these authors challenge the commonly held opinion 
that the strong coupling limit of QCD breaks chiral symmetry, irrespectively 
of Nf. In the already mentioned paper [9], the authors consider Nf = 16, the 
largest value of Nf still compatible with asymptotic freedom, and confirm that 
chiral symmetry is not broken at weak coupling by observing the expected bulk 
transition at intermediate coupling. Several other studies have observed phase 
transitions as a function of the gauge coupling for Nf = 8 [T5fl6jll7lll8jll9] . 
leaving the question of the fate of these transitions in the continuum limit 
open. Recent work by Appelquist and collaborators |20j finds evidence that 
Nf = 8 is still confining in the continuum limit. In particular, they conclude 



that Nf, the threshold of the conformal window, should be close to 12. They 
study the approach to the continuum by use of the Schrodinger functional and 
the step scaling method, the lattice discretized version of the renormalization 
group running coupling, and do not discuss how a bulk transition would affect 
these calculations. 

In summary, most of older studies at Nf = 8 suggested that the theory is 
already in the conformal window, while others confirmed analytic calculations 
[T|2|6f8] . indicating that Nf > 8, close to the upper limit for Nf calculated in 

i- 

The purpose of this study is to fully clarify the physics of eight flavours, as a 
first step towards the exploration of the phase diagram of QCD at large Nf. 



2 Simulation and Observables 



The simulations in this paper were performed using a slightly modified version 
of the publicly available MILC code [2T]. A setup similar to the one used by 
MILC in their recent paper on the QCD equation of state [22] was employed. 

We use an improved Kogut-Susskind fermion action, the "Asqtad" action 
which removes lattice artifacts up to 0(a 2 g 2 ) [2TJ and a one-loop Symanzik 
improved J22y23] and tadpole improved [21] gauge action. The complete action 
for Nf mass degenerate flavours can be written as 

Nf 
S = — — Tr In M(am, U, u ) + S gauge , (1) 



where 
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with M(am, U, uq) the fermion matrix for the Asqtad staggered action for a 
single flavour with mass m, and with a s = — 41ogwo/3.0684. The <Sj's contain 
all the lxl plaquettes, the 1x2 and 2x1 rectangles and the lxlxl 
parallelograms, respectively, that can be drawn on the lattice. The [/(C) 's are 
the traces of the ordered product of link variables along C, all divided by 
the number of colours. The tadpole parameter uq is defined in terms of the 
gauge invariant average plaquette as u = (U(C))\ C€S . Exploiting the rather 
low sensitivity of the plaquette value to finite volume effects, the u values 
for each run were determined by a self consistency procedure on 12 4 and 16 4 
lattices for every (3 investigated. 



It has been proposed recently by members of the HPQCD collaboration [25] 
to include the effect of dynamical quarks into the one-loop improvement of the 
gauge action, the required coefficients of which they determined to first order 
through lattice perturbation theory. While it is true that a complete 0(a 2 ) 
improvement should fully account for gauge and quark loop contributions, the 
reliability of a one-loop truncated perturbative contribution from dynamical 
quarks at large Nf becomes questionable. Since the coefficients found in [25] 
are sizable and linear in the flavour number, a resummation would probably 
be called for, especially in the case of large Nf. Barring this, there is a risk that 
including a truncated series would severely over compensate the actual quark 
loop effects and in fact worsen the 0(a 2 ) improvement, for values of Nf > 3. 
The scaling violations that are being observed for the heavy quark potential 
with the customary 'partially' improved action are not completely negligible, 
but well under control [26J . On the base of these considerations and the fact 
that observables employed in this study are less sensitive to scaling violation 
effects, we did not include these one-loop dynamical quark improvements. 

To generate configurations with eight degenerate dynamical flavours, we used 
the rational hybrid monte carlo (RHMC) algorithm [27], which allows for 
simulating an arbitrary number of flavours through varying the number of 
pseudofermions used, i.e. it inherently uses an analytical continuation in the 
number of flavours. Simulations were done with two pseudofermions, thus 
avoiding fractional powers of the Dirac operator in the molecular dynamics, 
and actually maintaining the best performance if compared with other choices 
of the number of pseudofermions. We are fully exploiting the characteristics 
of the RHMC algorithm and monitor possible effects of a rooting procedure - 
by the use of fractional powers of the Dirac operator - in ongoing studies at 
varying number of flavours [28J . 



To systematically investigate the character of eight flavour QCD under differ- 
ent thermodynamical conditions, simulations were run at two different values 
for the length of the temporal extent in lattice units N t , namely 6 and 12. As 
was mentioned in the introduction, the scaling behaviour of the critical lattice 
coupling, defined according to eq. (pj) as (3 C = 10/g 2 c , with this parameter is 



an indication of either the thermodynamical or bulk nature of the observed 
phase transition. To check for the influence of finite volume effects and deduce 
the scaling properties of observables, simulations were run for three spatial 
extents of the lattice N s = 12, 20, 24 for N t = 6, and N s = 24 for N t = 12. 
We chose a fixed value of the lattice degenerate quark mass am = 0.02 to 
explore the critical region, and will describe how results can be modified by 
doing simulations at a lower quark mass in section |4} As it will be clear from 
the outcomes, the chosen masses are sufficiently light to clearly disentangle 
the critical behaviour and establish the nature of the transition. 

The Monte Carlo history was collected with trajectories of total length r = 0.3 
to 0.4, and time step 5t = 0.003 to 0.007 from the lowest to the highest (3 
values and from the smallest to the largest volume. 



At each (3, the expectation value of the real part of the Polyakov loop 

1 . Nt 

L ~ ^v^ J^ReTr J] U 4 (x,x 4 ) 



(4) 



was determined. It is important to remember that the real part of the Polyakov 
loop is a true order parameter only for the pure gauge theory, which is recov- 
ered in the infinitely heavy mass limit m — > oo of the theory with dynamical 
flavours. It is not an order parameter in all other cases, though a clear change 
in its value and its susceptibility can be observed at sufficiently heavy or suf- 
ficiently light values of the fermion masses. Where the theory does enter the 
regime of 'light masses' is based on empirical observations [29], it might shed 
light on the mechanism relating confinement to chiral symmetry breaking, and 
it depends on the flavour content of the theory. 



The chiral condensate for Nf degenerate flavours in lattice units 
N f 
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was determined by using a stochastic estimator with 20 repetitions. The chiral 
susceptibility, measuring the variation of the chiral condensate with varying 
the fermion mass x = d{ifjifj)/dm at fixed /3 can be divided into a connected 
and disconnected component x — Xconn + Xdisc, given in lattice units by 
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respectively. We have conveniently written the condensate and its susceptibil- 
ities in terms of traces of (products of) the staggered fermion matrix M as 
they are actually computed in the simulation. The connected and disconnected 
contributions to the chiral susceptibility are measured separately; more on this 
can be found in [30]. The connected contribution can also be measured in a 
partially quenched manner, performing a numerical derivative of the chiral 
condensate with respect to the valence quark mass 

= d ML (7) 

Xconn dm v • [ ' 



Both measurement methods were implemented, with results being in excellent 
agreement. The complete set of values reported on here have been found using 
the first method, since it has a lower computational overhead. 

The disconnected chiral susceptibility is a non-local quantity that can be esti- 
mated from the variance of the bulk behaviour of the chiral condensate. How- 
ever, such a variance for a condensate computed with stochastic estimators 
will automatically include part of the connected contributions, through ran- 
dom sources multiplying themselves. Following Bernard et al. [30J , we straight- 
forwardly eliminate those contributions by only considering the off-diagonal 
elements of the covariance matrix of the random sources introduced for the 
estimation of the chiral condensate. 

We can use the chiral susceptibility and the chiral condensate to define two 
physically relevant quantities 

^ _ d{j}j}) _ 

Xu — X a Xconn T Xdisc \o) 



and 



m 



They are related through Ward identities [3T] to the spacetime volume integral 
of the scalar (er) and pseudoscalar (ir) propagators 



Xu.TT 



J d 4 x G^{x) , (10) 



thus implying that they should become degenerate when chiral symmetry is 
restored, following the degeneracy of the chiral partners. As a result, their 
associated cumulant R n = Xo-fXn should be, in the absence of explicit chiral 
symmetry breaking, one in a chirally symmetric regime, while it should be zero 



in the spontaneously broken phase. This dimensionless quantity is therefore a 
most useful physical observable in analyzing chiral phase transitions, and we 
will employ it to compare results from simulations run at different temporal 
extents of the lattice. 



3 Results at N t = 6 



The introductory discussion tells us that we must observe a phase transition 
between a phase where chiral symmetry is broken to a phase where chiral 
symmetry is realized, for some value of the lattice gauge coupling /3 C . There 
are basically two ways to assess a phase transition on a lattice. First, one 
can rely on obtaining infinite volume estimates of relevant quantities by per- 
forming simulations on lattices of different spatial extents. Whenever this is 
possible, the critical behaviour is analogous to the one observed in the con- 
tinuum. This strategy usually breaks down very close to the critical point, 
where the correlation length of the system grows large and becomes compa- 
rable to the lattice size. When this happens, it is mandatory to use finite size 
scaling techniques. The determination of the regime we are in can be made a 
posteriori, by comparing results obtained from different volumes. 

For the case at hand, results for the chiral condensate and the Polyakov loop 
are shown in Fig. [2j From the small differences found at different spatial 
extents, it was concluded that our results can indeed be considered infinite 
volume estimates for (3 < 4.1 and f3 > 4.15, and the /3 dependence is smooth. 
The jump between the two branches is very clear for both observables and 
suggestive of a discontinuity. 

These results allow us to immediately identify a critical region for j3 between 
4.1 and 4.15. In this region, one would expect to observe tunnelling between the 
two phases and enhancement of the associated susceptibilities. As mentioned 
above, the approach to the infinite volume limit should be slower here, because 
of the increasing correlation lengths. Indeed, volume dependence is observed, 
with the smaller of the lattices having a slightly higher (3 C . This behaviour 
is consistent with known predictions for finite volume effects on first order 
phase transitions (32]- Since the shift is much less pronounced for the change 
of volume between N s = 20 and 24, we concluded the infinite volume limit 
value for (3 C should be located at the lower parts of the designated critical 
region, and is bounded by 4.10 < (3 C < 4.125. For definiteness, we locate it at 
a value of 4.1125. 

The chiral susceptibilities are shown in Fig. [3] and confirm our picture. The 
residual splitting is associated with the explicit breaking induced by the mass 
term, but the tendency towards degeneracy is very clear. Finally, we plot the 



0.15-1 



0.10- 






0.05 



0.00 



og 



Q 



6 o A o g o 






o 



r 0.06 



-0.04 



-0.02 



0.00 



3.6 



3.8 



4.0 
P 



— I 1- 

4.2 






-0.02 



4.4 



Fig. 2. The chiral condensate (red) and Polyakov loop (blue) as a function of the 
lattice coupling (5. Measurements were done for three different spatial extents of the 
lattice N s : 12 (Q>), 20 (A) and 24 (O). The critical region has been indicated by 
vertical lines. 

cumulant R n , defined earlier, in Fig. [3} In the chiral limit, the cumulant should 
be zero in the broken phase, and approach one in the chirally symmetric phase. 
The observed trend is consistent with all other findings. 



3. 1 The critical region 



Having bracketed off the critical value for j3, we investigated the area in which 
the phase transition occurs in more detail. Of specific interest here is the 
potential presence of metastabilities. Simulations were performed on 12 3 x 6 
lattices starting from thermalized configurations at higher and lower values 
of /3, but with identical parameters. The presence of metastable states could 
show up, either through the presence of hysteresis effects, or by appearance 
of tunnelling effects in the Monte Carlo history. The latter effect would be 
expected to be somewhat enhanced, because of the usage of a smaller spatial 
extent. However, for all of the values of /3 that were examined near the thresh- 
old of the lattice of this size, a persistence of initial states was found after a 
large number of RHMC trajectories. This indicates the presence of a strong 
hysteresis loop with an extent equal to about half of our designated critical 
region, the consequences of which for the main observables have been plotted 
in Fig. |U 
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Fig. 3. The scalar Xa (blue) and pseudoscalar x-k (red) chiral susceptibilities as a 
function of the lattice coupling (3, confirming the degeneracy of the chiral partners 
in the symmetric phase. Measurements are shown for three different spatial extents 
of the lattice N s : 12 (Q), 20 (A) and 24(0).The cumulant R n is also shown (black, 
top of figure). The critical region has been indicated by vertical lines. 

Finally, it is a matter of interest to what extent the behaviour of the Polyakov 
loop, as a signal of confinement, is intertwined with that of the chiral con- 
densate. There are indications that the confinement/deconfinement transition 
is, in the presence of light quarks, driven by the chiral phase transition. Even 
though the signal of the Polyakov loop tends to be statistically weak for small 
numbers of measurements, it is possible to observe its change from the distri- 
bution of its phase. While the Polyakov loop tends to be distributed symmet- 
rically in the complex plane for confining configurations, whenever fermions 
are present in a deconfined system, it will tend towards the real phase of its 
Z(3) centre. In Fig. [5j we exhibit the Monte Carlo history of the phase of 
the Polyakov loop for a run at the threshold of the critical region and at the 
largest spatial extent, i.e. (3 = 4.1 and N s = 24. For comparison, the plot 
includes the equivalent part of the history at values of (3 before and after the 
transition. Note the presence of a metastability close to the threshold, also for 
this largest volume. This shows that the Polyakov loop, though it is a true 
order parameter only in the pure gauge theory, retains critical behaviour close 
to the transition, even in the presence of eight degenerate flavours. 

All data at N t = 6 point towards a first order phase transition. The location 
of the jump of the chiral condensate and the Polyakov loop at all three spatial 
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Fig. 4. Hysteresis effect for the chiral condensate (red) and Polyakov loop (blue) 
within the critical region, data are for the spatial extent N s = 12. Arrows indicate 
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broadest critical region which covers the jump of the chiral condensate and Polyakov 
loop at all spatial volumes and the hysteresis cycle. 
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Fig. 5. Monte Carlo history for (3 = 4.1 and N s = 24, at the threshold. The distri- 
bution flattens around zero in the deconfined phase, while it is randomly distributed 
between ±f in the confined phase. The same histories are shown for runs below and 
above the critical region (at (3 = 4.0 and 4.2) in lighter colours. The points are enu- 
merated according to their trajectory number, disregarding the first 100 trajectories 
after the cold start for each. 
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volumes considered, and the location of the hysteresis cycle at the smallest 
spatial extent, allows us to give the estimate of the broadest critical region, 
to be 4.10 < (3 C < 4.15. We take the occurrence of the jump of the chiral 
condensate and the Polyakov loop at our largest spatial volume N s = 24, 
as the best indicator of the critical point, giving f3 c (N t = Q,N S = oo) = 
4.1125 ± 0.0125, corresponding to the upper and lower bound 4.10 < /3 c (N t = 
6,iV s = oo) < 4.125. 



4 Results at N t = 12 



As explained in the introduction, the scaling of j3 c with varying the lattice 
temporal extension N t should distinguish between a thermal and a bulk tran- 
sition. A thermal transition occurs at a physical temperature T c which we 
determine as 

Tc = WWt (U) 



for a given temporal extent N t and critical lattice spacing a(/3 c ). Thus, per- 
forming simulations at different N t , the critical coupling should rescale in order 
to give the same critical temperature T c . Uncertainties in the determination of 
the critical temperature can be mainly introduced by lattice violations to the 
asymptotic scaling of observables. It is thus a possibility that past simulations 
leading to apparent scaling violations and indicating the presence of a bulk 
transition, might have been performed in a region of the lattice parameter 
space not sufficiently close to the continuum limit. Our highly improved lat- 
tice action should help us approaching the continuum limit and the correct 
scaling behaviour. 

We have thus performed simulations at N t = 12 with spatial volume N s = 24 
and at the same lattice mass am = 0.02, and looked for signals of the transition 
in the behaviour of the chiral condensate. The main result is shown in Fig. |6j 
where we observe a jump in the chiral condensate suggestive of a discontinuity, 
and slightly distorted and smoothed by the not so large spatial extent N s = 24 
as compared with the temporal extent N t = 12. Notice also that effects of 
explicit chiral symmetry breaking due to the non zero physical mass are to 
be more pronounced at N t = 12, for the same value of the lattice mass. 
The derivative of the condensate is also plotted in Fig. [6] to better locate the 
transition. The best fit value for the maximum of the derivative indicates a 
value for the critical coupling j3 c = 4.34±0.04. This result allows us to place an 
upper and lower bound on the critical coupling (3: 4.30 < (3 C < 4.38 at N t = 12, 
with iV s = 24. While finite spatial volume effects will not effect our conclusion, 
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Fig. 6. The chiral condensate at Nt = 12 and N s = 24 in lattice units as a function of 
the lattice coupling (3. Best fit curves are superimposed and the vertical lines indicate 
the critical region. The absolute value of the finite difference between measured values 
of the condensate, as an approximation of its first derivative, is plotted in the bottom 
part of the figure with an arbitrary rescaling. It shows a peak at (3 = 4.34. 

we cannot quote this value as the infinite volume value, and we postpone a 
more refined analysis of finite volume effects to a future publication. 

The combined results at N t = 12 and N t = 6 strongly suggest the occurrence 
of a first order thermal transition. To confirm this, an asymptotic scaling 
analysis is needed in order to verify that we are actually measuring a critical 
temperature in the continuum real world. We do so by means of the standard 
relation that connects the lattice cutoff A^ to the gauge coupling g, aA^ = 
R(g 2 ) with 



Rtf) = (b g 2 



-bl/2&2 p -l/2b g 2 



(12) 



where the two loop RG running of the /3-function is accounted for, with the 
universal one- and two-loop coefficients given by 



bo 



1 



16tt 2 
1 

(16tt 2 ) 2 



11 



3' 
102 



-Ni 



38 



N< 



(13) 



14 



for Nf massless flavours. Given a physical temperature T c , eq. (11) implies the 
scaling relation N t R(g c (N t )) = const. Solving for N t = 6 and N t = 12, we can 
predict (3 c (g c (N t = 6)) by knowing (3 c {g c {N t = 12)). A strong discrepancy with 
the actual lattice determination might besuggestive of a bulk zero temperature 
transition, while a small discrepancy can indeed be expected and imputed to 
violations of asymptotic scaling and residual effects due to a non zero fermion 
mass. The question is what coupling /3 c (g c ) is the most appropriate one to in- 



sert in eq. ( 12 ) in order to compare continuum perturbation theory with lattice 
calculations. Our simulations use a one-loop Symanzik improved and tadpole 
improved gauge action. These improvements are expected to ameliorate the 
agreement with the asymptotic scaling by correcting for up to and including 
0(a 2 ) effects. In a finite temperature study this translates into the statement 
that the agreement with asymptotic scaling will set in at lower values of N t 
with respect to an unimproved lattice action. 

Given these premises, we have predicted the critical coupling (3 c {N t = 6) by 
using an effective coupling which relates the lattice (3 and g 2 as (3 = 6/g 2 



in eq. (12). Notice that the lattice bare coupling, coefficient of the plaquette 
action in eq. (IT]), is (3 p i = 10/ g$ where the rescaling of the coefficients is due to 
the one- loop Symanzik improvement. The use of (3 — Q/g 2 in the RG formula 
is meant to effectively take into account the tadpole and Symanzik one-loop 
improvement in our action. 

By using (3 c (N f = 12) = 4.34 ± 0.04 we obtain the prediction (3 c (N t = 6) = 
4.04 ±0.04, which deviates by less than 2% from the lattice determined f3 c . We 
also verified that a rescaling of the effective coupling (3 — > (3uq 4 [21] improves 
the prediction by only 0.5% indicating that our effective coupling prescription 
is in fact already accounting for all the improvement. Further improvement of 
the lattice gauge action - e.g. the proper inclusion of quark loop corrections 
in the improvement program - and corrections for the non zero fermion mass, 
should account for the tiny residual asymptotic scaling violations in this study. 

In Fig. [7] we compare the transition at N t = 6 and at N t = 12 by means of the 
first derivative of the chiral condensate as a function of T/T c as predicted by 
the asymptotic scaling using the effective coupling (3 = Q/g 2 - As before, the 
lattice determined l3 c (N t = 12) corresponds to the point T/T c = 1. The most 
important observation is that the transition at N t = 6 is close to the predicted 
point T/T c = 1, actually happening at T/T c ~ 1.20. We show the data of the 
first derivative with a gaussian fit superimposed. For the data at N t = 6 
we show, for the sake of comparison, the curves at the largest spatial extent 
N s = 24 and the smallest N s = 12. The lattice determined critical regions at 
Nt = 6 and Nt = 12 overlap at their boundaries. The amount of non-overlap 
is a measurement of the residual but small violations of asymptotic scaling. 
We did not observe any signal of discontinuity in the Polyakov loop in our 
simulations at N t = 12, with a lattice mass am = 0.02 and in the critical 
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Fig. 7. 'Scaling plot' for the finite difference approximation to the absolute value of 
the first derivative of the chiral condensate as a function ofT/T c , determined using 
the effective coupling (3 = 6/g 2 in the RG formula. Data at Nt = 6, N s = 24, 12 
(blue) are compared with data at N t = 12, N s = 24 (red). The point T/T c = 1 cor- 
responds to (3 c {Nt = 12) = 4.34. A gaussian fit is superimposed to the Nt = 12 data, 
while indicative gaussian curves are shown for Nt = 6. The baseline is subtracted 
and we only indicate data near the critical points. The lattice determined critical 
regions at Nt = 6 (blue) and Nt = 12 (red) are indicated by vertical lines. 



region of the chiral condensate. While the presence or absence of this signal is 
of no physical relevance below the infinite mass limit, given that the Polyakov 
loop is not an order parameter in the presence of dynamical fermions, we 
can still use its behaviour at light quark masses to gain further insight into 
the mechanism that relates chiral symmetry restoration to deconfinement and 
the two limiting regimes of zero mass and infinite mass fermions. Here, we 
limited ourselves to verify that the presence or absence of a discontinuity of the 
Polyakov loop in coincidence with the discontinuity of the chiral condensate 
depends upon the physical fermion mass. By sufficiently lowering the fermion 
mass at Nt = 12 we expect that the picture of Nt = 6 will be restored. 

However, the region of sufficiently light quark mass is a dynamical property 
of the Nf = 8 theory which cannot be judged a priori. A further refinement 
of this study can help in locating the critical end-points. The reader might 
want to review |29j for a pioneering study at Nf = 4. We postpone to future 
publications a more refined study of the quark mass dependence of the chiral 
condensate, and a thorough scale setting procedure for the results obtained in 
this work. 
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The present results allow us to firmly conclude that we are seeing a true 
thermal transition, in other words Nf = 8 undergoes a chiral restoration 
transition at finite temperature. 



5 Discussion 



We find (3 c {N f = 8, N t = 6,N S = oo) = 4.1125 ± 0.0125, or 4.1 > (3 c {N f = 
8, N t = 6, N s = oo) < 4.125 and (3 c {N f = 8,N t = 12, N s = 24) = 4.34 ± 0.04, 
or 4.30 > p c (N f = 8,N t = 12, N s = 24) < 4.38 and observe the asymptotic 
scaling of the critical temperature within 20% by use of the perturbative scal- 
ing induced by asymptotic freedom, where the two loop beta function for an 
SU (3) gauge theory with Nf = 8 massless flavours has been used. These ob- 
servations, obtained with an improved lattice action, might well explain the 
lack of asymptotic scaling noted in previous studies, where no improvement 
was yet considered. 

All evidence presented in this paper is consistent with asymptotic scaling and 
the occurrence of a true thermal transition. The emerging picture confirms the 
conclusion of the recent work [2D] , as well as ladder calculations [8], providing 
a strong evidence that SU(3) gauge theory with eight flavours is in the normal, 
chirally broken phase of QCD at zero temperature and in the continuum limit. 
This study is one step in our attempt to clarify the way SU(3) gauge theory 
with many flavours approaches the conformal phase [28] that precedes the loss 
of asymptotic freedom. 
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